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Abstract 

A spherical quadrangulation is a loopless graph embedded on the 
sphere such that each face is bounded by a walk of length 4, parallel 
edges allowed. The family of the isomorphism classes of quadrangulations 
can be generated by a sequence of graph operations called vertex splitting, 
starting from the path graph P2 with three vertices and two edges. Pi is 
also referred to as the common ancestor of all quadrangulations. We define 
the degree 1 < D < |_cZ/2j of a splitting 5 (where d is the degree of the split 
vertex) and consider restricted splittings Sij with l<i<D<j< [d/2\ . 
As Brinkmann and coworkers have recently pointed out, restricted split- 
tings 52,3 generate all simple quadrangulations. 

Here we investigate the cases Si, 2, 5*1,3, Si, 1, 5*2,2, 53,3. First we show 
that the restricted splittings Si, 2 are exactly the monotone ones in the 
sense that the resulting graph contains the original as a subgraph. Then 
we proceed to show that they define a set of nontrivial ancestors beyond 
P2 and each quadrangulation has a unique ancestor. 

Our results have a direct geometric interpretation in the context of 
mechanical equilibria of convex bodies. The latter can be defined as the 
scalar distance R(6, if) measured from the center of gravity and the Morse- 
Smale complex associated with the gradient of R corresponds to a 2- 
coloured quadrangulation with independent set sizes s, u. The numbers 
s, u of coloured vertices identify the primary equilibrium class associated 
with the body by Varkonyi and Domokos. We show that the 5i,i and 52,2 
splittings generate all primary equilibrium classes (in case of 5i,i from a 
single ancestor, in case of 52,2 from a finite nontrivial set of ancestors). 
This is closely related to the geometric results of Varkonyi and Domokos 
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where they show that specific geometric transformations can generate all 
equilibrium classes. 

If, beyond the numbers s,u, the full topology of the quadrangulation 
is considered, we arrive at the more refined secondary equilibrium classes. 
As Domokos, Langi and Szabo showed recently, one can create the ge- 
ometric counterparts of unrestricted splittings to generate all secondary 
classes. Our results show that restricted, monotone splittings Si, 2, while 
adequate to generate all primary classes from one single ancestor, can 
only generate a limited range of secondary equilibrium classes from the 
same ancestor. The geometric interpretation of the additional ancestors 
defined by monotone splittings shows that minimal polyhedra play a key 
role in this process. We also present some computational results on the 
cardinality of secondary equilibrium classes and multiquadrangulations. 

Keywords: plane multiquadrangulation ; unrooted unsensed colored map ; 
vertex splitting ; census ; convex body ; equilibrium class 



1 Introduction 

Our work is motivated by the classification system of convex, homogeneous 
3D bodies introduced by Varkonyi and Domokos [20 . They map each body 
its primary equilibrium class defined by the numbers of the stable and unsta- 
ble equilibrium points of the body surface. Moreover, the isomorphism classes 
of the topologies of the equilibria provides a refined, secondary classification 
system, where such a topology can be genuinely represented by a 2-colourcd 
quadrangulation [7]. 



1.1 Generating multiquadrangulations 

A quadrangulation of the sphere is a loopless graph embedded in the sphere 
having every face bounded by a closed walk of length 4. We allow parallel 
edges, and the boundary walk may repeat edges or vertices. This definition was 
also used by Mohar et al. [TS], however, Archdeacon et al. [T] applied the word 
"pseudoquadrangulation" instead for multigraphs. If we want to emphasize that 
the quadrangulation may have parallel edges, it is called a multiquadrangulation, 
if a quadrangulation has no parallel edges, it is called a simple quadrangulation. 
Note that the 2-path P2 (the path of length 2 with two edges and three vertices) 
is the smallest quadrangulation, and the 4-cycle C4 (the cycle of length 4) is the 



smallest simple quadrangulation, illustrated on Figure 5| 

Two quadrangulations are considered isomorphic, if there is a homeomor- 
phism from one to another that cither preserves or reverses the orientation of 
the embedding. Equivalently, the cyclic ordering of the incident edges at each 
vertex is either preserved or reversed, so e.g. a graph is isomorphic to its reflec- 
tion. Such an isomorphism class is also called an unsensed, unrooted map in the 
literature [22]. Let J3 denote the family of all multiquadrangulations, and £H\ 
the family of all simple quadrangulations. 
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Figure 1: Vertex splitting of degree m. 



We say a graph family & is generated from the starting set K C & by some 
given graph operations, if each graph in & can be constructed from some graph 
of K by applying a finite series of the given graph operations. All our graph 
operations are based on the vertex splitting, depicted on | Figure T\ explained 
as follows. The embedding (the cyclic ordering) of the graphs showed is also 
important so the small open triangles denote that other edges may occur only at 
that position so they clarify the cyclic ordering of the edges in the original (left 
to the arrow) and in the resulting (right to the arrow) graph. Vertex splitting 
is also known in the literature as the inverse of face contraction. 

Batagelj ,3], Negami and Nakamoto [15] showed that the splittings generate 
£?i from C4. There are a number of related results regarding inductive genera- 
tion of certain simple quadrangulation families, e.g. Batagelj [5] gave inductive 
definition of 3-connected quadrangulations, Nakamoto |17j generated quadran- 
gulations with minimum degree 3, Brinkmann et al. [6] improved efficiency for 
all these families above. While [3J [TBI HH] mainly focus on simple quadrangula- 
tions, the following observation is a straightforward extension of their results: 



the splitting generalized for parallel case as well (m = 1 on Figure 1(b) I gener- 
ates =2 from P 2 . 

As shown on Figure 1 the splittings replace a vertex v with vertices w and 
v' dividing the edges of v. The degree of a splitting is D := mm{d(v'),d(w)} 
where d(v) denotes the degree of the vertex v. Note that 1 < D < \_d{v)/2\ 
always holds. We consider restricted splittings Sij with a limited range of D 
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Figure 2: The monotone vertex splittings Si,2- The 2-splitting has two variants: 
simple case and parallel case. 



such that l<i<D<j< [d/2\. Brinkmann et al. [5] showed that the 
restricted splittings 6*2,3 (with 2 < D < 3) are enough to generate £l\ from 
C4. Their result can be easily extended for parallel graphs using the generalized 
vertex splitting, so our previous observation can be improved as follows: 

Theorem 1. The restricted splittings S*i,3 generate J2 from P^. 

In this paper we investigate the hierarchy generated by restricting the split- 



tings Si t 2 (see Figure 2 ), and we also show some application of restrictions 61,1, 
52,2 a n d S3, 3- The main reason we focus on splittings S*i,2 is that these are not 
only local modifications but they purely extend the graph without removing any 
edge. Formally, we say a graph operation is monotone if the original graph is 
the embedded subgraph of the resulting one. Actually Si 2 are the only mono- 
tone operations on quadrangulations introducing only one new vertex. The local 
monotone operations are also called face subdivisions in the literature meaning 
one face is divided into smaller regions (e.g. (HE])- Tutte [T§] called the order 
of a face subdivision the number of the introduced vertices in the operation, 
although he applied this concept on triangulations. So on quadrangulations the 
monotone vertex splittings are exactly the face subdivisions of order 1. 

For the splitting Sd.d with a given D, we use the shorthand D-splitting, and 
we call its inverse D- contraction. If a D-contraction is applicable to a graph, we 
say the graph is D-contractible, otherwise we say the graph is D -irreducible. It 
is known that among the simple quadrangulations, C4 is the only graph which 



is D- irreducible for any D [3J [TS] . It follows from Theorem 1 that among the 
multiquadrangulations, only P2 is D-irreducible for any D, hence P2 is also 
referred to as the only ancestor of all quadrangulations. However, restricting 
the splittings to 5*1,2 also admits other, nontrivial ancestors. Throughout this 
paper we say a graph is an irreducible ancestor (or shortly irreducible) if it 
is 1-irreducible and 2-irreducible. We mention that the concept of irreducible 
graphs with respect to some given graph operations was analogously used in e.g. 
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pa m. 

We characterise the irreducible graphs as the ones with minimum degree 3, 
and extend a result of Batagelj [3] that the small graphs with less than 8 vertices 
are not irreducible, except P 2 : 

Theorem 2. Every quadrangulation with less than 8 vertices is generated from 
P2 by monotone splittings. 

We strengthen this theoretical result with observations on a data set gen- 
erated by our computer program which we developed to explicitly enumerate 
every possible quadrangulation for a given size, based on the software plantri 
j5]. Furthermore, our data set also shows that there are only three irreducible 
graphs up to a size of 10 vertices. Besides P2, the other ones are the radial 
graphs (a.k.a. vertex-face incidence graphs [T3]) of the skeletons of the two 
smallest polyhedra. The radial graph R(G) of an embedded graph G is a bi- 
partite embedded graph such that one of the independent vertex sets of R(G) 
corresponds to the vertex set of G, the other one to the face set of G. Two 
vertices are connected in R(G) with the same multiplicity as the incidence mul- 
tiplicity of their preimages in G (i.e. the appearance count of a vertex in the 
boundary walk of a face), with the obvious cyclic orderings. In addition, we 
prove 

Theorem 3. The radial graph of any polyhedral skeleton is irreducible, and they 
are generated from the radial graphs of the pyramids with polygonal base by the 
restricted splitting 83,3. 

Then we prove that the irreducible ancestor of any graph is unique, i.e. a 
graph cannot be generated starting from two different irreducible graphs: 

Theorem 4. If a quadrangulation can be generated from both irreducible an- 
cestors Ai and A 2 by monotone splittings, then Ai is isomorphic to A 2 . 



1.2 Generating secondary equilibrium classes 

The main goal of this paper is twofold: one is to contribute some observations 
on the reachability of quadrangulations i.e. which quadrangulations can be gen- 
erated from others, the other is to apply them on the geometric interpretation 
mentioned earlier. So finally, we translate our observations into the context of 
topology of the equilibrium points of convex bodies. A generic convex body is 
given by its scalar height function R(0, ip) which gives the distance between the 
surface and the mass center for any direction. A surface point is an equilibrium 
if the gradient of the height function R is zero at that point (illustrated on the 



ellipsoid on Figure 3 ) . An equilibrium is stable if R takes a local minimum 
value, unstable if R takes a local maximum value, otherwise it is a saddle point. 
According to the Poincare-Hopf theorem [2J, the cardinality of stable, unstable 
and saddle equilibria of a body, denoted respectively by s, u, h, are deeply re- 
lated: s + u — h = 2, thus the pair {s, u} defines unambiguously the primary 
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(a) Height function R 



(b) Gradient field of -R 




(c) Heteroclinic orbits 





(d) Height function R 



(e) Gradient field of —R 




(f) Heteroclinic orbits 



Figure 3: Morse-Smale graph of the ellipsoid. The axes of the ellipsoid are 
rotated for presentation purposes. Colouring notation: red: stable, blue: un- 
stable, white: saddle. 



6 



equilibrium class of the body [50] • The primary equilibrium classes can be re- 
fined into secondary equilibrium classes which determine also the topology of 
the equilibria defined by the Morse Smalc complex of the height function R [7] . 
This topology is usually given by a 3-coloured spherical quadrangulation, where 
the vertices are the equilibria, the colour of a vertex defines its equilibrium type 
(stable, unstable or saddle) , and an edge is a particular path on the surface con- 
necting two equilibria, whose tangent vectors agree with the gradient vectors of 
R ( |Figurc 3 ). We transform these 3-coloured surface graphs with a bijection to 



generic 2-coloured quadrangulations keeping the underlying geometric meaning, 
and we use this latter form to define the secondary class as a fully combinatorial 
object. We mention that a spherical quadrangulation is always bipartite [Tl [T2"] 
hence the 2-colouring is possible. 

The data set yielded by our program enables us to present some statistics on 
number of possible secondary classes as well. Because of the 1-1 correspondence 
of 2-coloured quadrangulations and connected plane graphs [TTJ H>] , our num- 
bers agree with the numbers of the unrooted and unsensed maps calculated by 
Wormald [H] , and exhaustively enumerated by Walsh [2TJ H3] • Walsh [52] sur- 
veys different types of map census results achieved both exhaustive search and 
with formulae. We are also able to give statistics on the number of multiquadran- 
gulations using a trivial relation between the number of maps, quadrangulations 
and self-dual quadrangulations for a fixed size. 

Varkonyi and Domokos |20| constructed the geometry of a representative 
body of the primary class {1,1} referred as mono-monostatic body, also known 
as Gomboc. They also concluded that every primary class for all s,u > 1 is 
generated from the Gomboc via their specific geometric transformations called 
Columbus' algorithm. Columbus' algorithm is a sequence of sensitive modifica- 
tions on a body perturbing the surface only at the vicinity of an equilibrium, 
such that the body belonging to the primary class {s, u} is transformed to an- 



other one belonging to {s + l,u} or {s,u + l} (illustrated on Figure 4). Because 
of the geometric feasibility of Columbus' algorithm, i.e. the transformations al- 
ways can be applied around any equilibrium, they referred to the Gomboc as 
the ancestor of every primary class. 

The combinatorial equivalent of the original Columbus' algorithm is the 
sequence of monotone coloured splittings. The result of [5D] can be reformulated 
as follows: for any s, u, some secondary class in the primary class {s, u} is 
generated from the Gomboc using the monotone coloured splittings. (We detail 



the secondary class of Gomboc in section 4] ) We extend this statement with the 



following corollaries of our theorems: 

Corollary 1. Not every secondary class is generated from the Gomboc using 
the monotone coloured splittings, but the ones with s + u < 8 are all generated 
from the Gomboc. 

We say a polyhedron is a minimal polyhedron if its every face contains one 
stable and its every vertex is an unstable equilibrium. With respect to the 
monotone coloured splittings (and the original Columbus' algorithm), they are 
also ancestors: 
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(a) Cutting the surface with a plane near Uq 




(b) Modified gradient field of — R 



Figure 4: Applying a step of Columbus' algorithm on an egg-shaped surface. 
Two equilibria appear on the surface: Si , Hi . 

Corollary 2. The secondary classes of the minimal polyhedra are irreducible 
so not generated from any other secondary class. 

Finally, if we aim to generate a secondary class from every primary class, it 
is enough to use one of the two monotone coloured splittings: 

Corollary 3. It is enough to use either the coloured splitting S± t i or to 
generate every primary class from a finite starting set of secondary classes. 

We mention that the coloured splittings are also a restricted subset of the 
operations called "cancellations" in computational geometry by Edelsbrunner 
et al. |10) and Bremer [4, Fig. 7.3], who used them to simplify a multi-resolution 
mesh structure. Domokos, Langi and Szabo [7] showed recently that the geo- 
metric counterparts of the unrestricted coloured splittings generate the whole 
family of secondary classes, however, this is beyond the scope of this paper. 

1.3 Organization of this paper 

The paper is organized as follows. Section 2 analyzes the vertex splitting and 
establishes its properties on the reachability of quadrangulations. Section 3 
presents the hierarchy of irreducible quadrangulations: the reachability of small 
graphs, the polyhedral irreducible quadrangulations and the uniqueness of an- 
cestors. Section 4 interprets our results on the equilibrium topologies of convex 
bodies. Finally we show some statistics on our data set in section 5. 
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l n -0 

(a) P 2 (b) C 4 (c) Q 3 (d) Qa 

Figure 5: The four smallest quadrangulations with n < 4. 

2 The properties of vertex splitting 

In order to generalize the splitting for parallel case (with the notations illustrated 



on Figure 1), let cr(v) denote the cyclic ordering of the vertex v in clock-wise 
order, so a(v) = (ei, . . . , es v \\. Let rii denote the other endpoint of the edge 
e». 

Definition 1. A vertex splitting on a quadrangulation G is a graph operation 
transforming G to G', specified by a walk n\e\ve m n m of G. The vertex v of G is 
replaced by vertices w and u' in G' dividing the edges of v. The cyclic ordering 
of w is a(w) = (e' 1; . . . , e' m ), where e\ has the same other endpoint as had; 
the cyclic ordering of v' is <r(v') — (ei, e" l7 e m +i, . . . , e^)), keeping some edges 
of the former vertex v, and ej^ = e m if m > 1, otherwise ej^ has the same other 
endpoint as e' m and e m had. 

Note that if m > 1 and m is identical to n m then e\ and e m are distinct 
parallel edges. The new face introduced by the splitting is bounded by the 
walk v' einie^we^rime'^v' . Observe that the degree of the new vertices are 
d(w) — to and d(v') — d(v) — m + 2. In addition, the splitting specified by 
the walk n\e\ve m n rn creates the very same graph that the one specified by the 
walk n m e m ve\n\. Hence the degree of a splitting, denoted by D, is invariant to 
reflection, and D := mhx{d(v ), d(w)} = min{m, d(v) — m + 2}. 

Let n denote the number of vertices of the quadrangulation. There are 
exactly four quadrangulations such that n < 4, shown on |Figure 5[ where C4 is 
generated from P2 with a 2-splitting, and Q3, Q4 are generated from Pi with a 
1-splitting. For larger graphs, we have 

Proposition 1. For a quadrangulation such that n > 4, 

(i) if it has a vertex of degree 1, then it is 1-contractible even for n = 4, 

(ii) if it has a vertex of degree 2, then it is 2-contractible, 

(Hi) if its minimum degree is k, then it is k-contractible, and not l-contractible 
for any I < k. 

(iv) it is irreducible if and only if its minimum degree is 3, 
Proof. Part (i) is a straightforward consequence of the quadrangulated faces. 
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Figure 6: 3-splitting. The half edge at v denotes that there must be an edge 
otherwise it would be a 2-splitting. 

For part (ii), suppose a vertex w is incident to two edges going to n\ and ri2- 
If n\ and 712 are distinct vertices, then the graph is 2-contractible as depicted 
on |Figurc 2| If m = and n\ was not incident to any other parallel edge, the 



graph would be isomorphic to Q3 (see Figure 5 1 . If there are other parallel edges 



incident to n\, then the graph is 2-contractible, the parallel case is applicable 



(see Figure 2 1 



Part (iii) follows from the definition of the degree of splitting and from the 
observation that a fc-splitting introduces a vertex of degree k. 

Part (iv) follows from (i)-(iii) as the monotone splittings introduce a vertex 
of degree 1 or 2. □ 

[Proposition l| provides the necessary observations to prove that the splittings 
Sx t $ generate J2. 



Proof of Theorem 1. The statements follows from Proposition 1 and from the 
fact that the minimum degree of a multiquadrangulation is either 1, 2 or 3 by 
Eulcr's formula. □ 

It is clear that the 1-splitting and the 2-splitting are monotone, and modify 
only one face of the graph in the vicinity of a point. Hence the perturbation 
of the graph is minimal. Monotonia ty is also sufficient to characterize these 
splittings: 

Proposition 2. Every plane graph operation defined on quadrangulations which 
adds one vertex and is monotone, is a restricted splitting S\ t 2- 

Proof. Consider any operation satisfying the conditions, putting the new vertex 
w into an original quadrilateral face bounded by the walk X1X2X3X4.X1, listing 
only the vertices of the walk. Note that some of the vertices Xi may coincide. 
By Euler's formula, we need to add two edges in addition to the new vertex, so 
connect, without loss of generality, w and x\ with an edge. Then, we have an 
"almost quadrangulation" : except for one face which is bounded by the walk 
X\WX\X2XzXiX\ of length 6. This face can be divided into two quadrilateral faces 
by adding one edge in three ways: connecting x\ and X2 (the fourth and the 
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(a) Radial graph (left) and the skeleton 
(right) of the tetrahedron 




(b) Radial graph (left) and the skeleton (right) 
of the square pyramid 



Figure 7: Irreducible quadrangulations. 

seventh elements of the walk sequence), or x% and X4 (third and sixth elements), 
or w and £3 (second and fifth) with an edge. The first and the second case is 
applying a 1-splitting, the third case is applying a 2-splitting. □ 



3 Irreducible quadrangulations 

As being irreducible for n > 3 is equivalent to having minimum degree of 3, one 
can easily verify that e.g. the radial graph of the skeleton of the tetrahedron 
is irreducible (see Figure 7(a) [ )■ According to Theorem 2 it is the smallest 



irreducible quadrangulation after Pi. In preparation of proving it, we need to 
generalize slightly an earlier result: 

Proposition 3 (Batagelj [3]). Every simple 3-connected quadrangulation has 
at least 8 vertices of degree 3. 

The generalization of Batagelj 's proof for our need is straightforward as his 
statement still holds for multigraphs and instead of 3-connectivity it is enough 
to assume the minimum degree is 3. The latter observation has already been 
made in [B]. Hence we use the following statement: 

Proposition 4. Every quadrangulation with minimum degree 3 has at least 8 
vertices of degree 3. 

Proof. The layout of this proof is similar to Batagelj 's proof, but the conditions 
are relaxed. From Euler's formula we have that the sum of the degrees of a 
quadrangulation of size n := s + u is An — 8. If it has k vertices of degree 3 and 
n — k other vertices of degree at least 4, then the sum of the degrees is also at 
least 3k + 4(n — k), implying k > 8. □ 
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Proof of Theorem 2. By |Proposition 1 and Proposition 4 we have that the only 



irreducible ancestor such that n < 8 is P 2 . So the theorem can be proved by 
induction. □ 



The data set generated by our program has strengthened Theorem 2[ as we 



could count easily the irreducible ones for n < 10 by [Proposition 1 The statis- 



tics showed that there are exactly three irreducible graphs among them: P 2 and 
the radial graphs of the two polyhedral skeletons shown on |Figure 7| |Theorem 3| 
states that the radial graphs of polyhedral skeletons are all irreducible. More- 
over, it gives an inductive definition for them, as they cannot be generated by 
monotone splittings. 

Proof of Theorem 3. Any polyhedral skeleton is simple and 3-connected by Steinitz' 
theorem [T3]. Let i?4 denote the family of simple and 3-connected quadrangu- 
lations having no separating 4-cycles. It is known [3] that a graph is simple and 
3-connected if and only if its radial graph belongs to i?4, proving the first part 
of the theorem. 

It was also shown [6] that is generated by the restricted splitting £3,3 
starting from the pseudo-double wheels. A pseudo-double wheel is a cycle of 
even length, with its inner and outer face subdivided by a vertex, such that the 
inner vertex is adjacent to the odd-numbered vertices of the cycle, the outer 
vertex is adjacent to the even- numbered vertices of the cycle [5]. It is easy to 
verify that a pseudo-double wheel of size n = 2k + 2 is the radial graph of the 
skeleton of a pyramid with a fc-sided base which completes the proof. □ 

Now we show that the irreducible ancestor of a quadrangulation is unam- 
biguously defined, so a quadrangulation cannot be generated from different ir- 
reducible ancestors at the same time. Consequently, for any two different irre- 
ducible quadrangulations A\ 7^ A 2 the families generated starting from {Ai} 
and from {A2} are disjoint, thus the set of irreducible quadrangulations leads 
to a natural partition of J2. 

Proof of Theorem 4- Without loss of generality, suppose Ai is not isomorphic 
to P 2 thus the minimum degree of A\ is 3. Let Go, G\, . . . , G p be the series 
of graphs yielded by the generation process of G from A% where Go = A\ and 
G p = G, and, indexed in reverse order, G p + q , . . . , G p +i, G p the graphs yielded 
by the generation process from A 2 where G p+q — A 2 , for some p, q > 0. We 
will prove that the edges of Go remain intact despite applying any monotone 
vertex splitting or monotone face contraction, i.e. Go is the embedded subgraph 
of every G& for all < k < p + q, so A\ is the subgraph of A 2 , and vice versa. 

We say an edge of a graph Gi is ancient if it is also in Go- By induction, 
assume that Go is the embedded subgraph of the graphs Go, ... , Gfc_i for some 
k, which trivially holds for k = 1. If Gfe is created from Gfc_i by a monotone 
splitting, clearly no edge is removed due to the definition of monotonicity. If 
Gfc is created from G^_i by a 2-contraction, then we use a basic property of the 
2-contraction that if it removes an edge, then one endpoint of the edge is also 
removed and the degree of the removed endpoint is 2. However, according to 
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the induction hypothesis, each ancient edge of Gk-i has endpoints with degree 
at least 3 thus no 2-contraction is applicable to Gk-i to remove an ancient edge. 

The last case is where Gk is created from Gk-i by a 1-contraction, removing 
two edges: e[ and e" (see the definition of vertex splitting and Figure 1(b)). 
Edge e'j cannot be ancient because of having and endpoint of degree 1 as ex- 
plained in the previous case. If edge e'{ is ancient and ei is not, then again no 
ancient edge is removed as it is only a technicality if we actually remove the 
edge e'l or e\ because of the symmetry of the 1-contraction. If both edges e'{ 
and ei would be ancient and e[ would not be, then, because of the induction 
hypothesis, there would be a face in Go bounded by the walk v' e\n\e'[v' con- 
tradicting the assumption that Go is a quadrangulation. Thus a 1-contraction 
cannot remove an ancient edge either. 

Consequently, as neither the splittings nor the contractions alter the embed- 
ding of the rest of the graph, Go remains the embedded subgraph of Gk for all 
k. □ 



4 Generating secondary equilibrium classes 

In this section we interpret our results on generating secondary equilibrium 
classes as coloured surface graphs of convex bodies. Some of these observations 
were already outlined in the conference version of this paper [IB] . 

The idea of representing a 3D body with some surface graph appears in mul- 
tiple disciplines. Most famous are the polyhedral graphs that are the skeletons of 
the convex polyhedra, characterized by Steinitz's theorem [14] . Another related 
appearance is the mesh generation of physical shapes in computational geome- 
try, when some Morse-Smale complex of a body is drawn as a surface graph, see 
e.g. [S]. However, the underlying function used for meshing is not necessarily 
the height function R, thus the nodes of the mesh and the equilibrium points 
do not necessarily coincide. 



4.1 Geometric interpretation of quadrangulations and ver- 
tex splittings 

To determine the secondary class of a body, we need to introduce some concepts 



from Morse theory [51 [TU] , illustrated on Figure 3 In generic case, we say a path 
on the surface is a heteroclinic orbit, if its tangent vectors agree with the gradient 
vectors of the height function R, and its endpoints are two equilibria of different 
type (see Figure 3(c) and |3(f)[ ). It is known the heteroclinic orbits incident to 
a saddle point are isolated on the surface, and there are only a finite number of 
them. 

These isolated orbits divide the body surface into quadrilateral cells (see 
Figure 8(a)), and in each cell an infinite number of non- isolated, heteroclinic 



orbits are going from the unstable to the stable point which we disregard for now. 
In this way, the body surface defines a vertex-coloured multigraph embedded 
on the sphere, where the vertices are the equilibria, the edges are the isolated 
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(a) Morse-Smalc graph (b) Triangulation (c) Quasi-dual 

(isolated orbits) 

Figure 8: Determining the secondary class of the ellipsoid. 



heteroclinic orbits connecting saddle and non-saddle vertices, the faces are the 
quadrilateral cells, and the colour of a vertex gives its type of equilibrium. 

This graph is referred as the Morse-Smale graph, which means possessing 
three properties [2HTU]. (1) The graph is a quadrangulation of the plane unless 
the body is in the primary class of the Gomboc. (2) Every quadrilateral bound- 
ary walk is a sequence of a saddle, a stable, a saddle and an unstable vertex. 
(3) The degree of every saddle vertex is 4. 

We perform two invertible transformations on the Morse-Smale graph, to 
yield a more compact graph which is still a genuine representation of the geom- 
etry. First, we connect the stable and unstable points in each quadrangulated 



face, creating a triangulation (see Figure 8(b) ). The geometric interpretation of 



this operation could be that from each cell, we pick one orbit arbitrarily from 
the infinite set of unstable-stable orbits. Then we remove the saddle points and 
the edges incident to them, creating another quadrangulation leaving only the 



edges created in the previous step (see Figure 8(c) ). This idea also appears in 



[8] as the quasi-dual of the Morse-Smale complex, however, as their goal is to 
simplify of the mesh structure, they do not mention that it can be shown that 
these transformations are invertible, so we lost no information by removing the 
saddles. We use the resulting quasi-dual graph to define the secondary class of 
the body, summarized in the following 

Definition 2. A quasi- dual is a 2-coloured multiquadrangulation, where the 
independent sets are called the stable and the unstable equilibria, of size s 
and u, respectively, if s + u > 2. A quasi-dual of a body is the quasi-dual 
obtained with the process described above. For technicality, in case s + u = 2, 
we define the quasi-dual of the Gomboc as the path Pi of length 1 (with one 



edge) connecting a stable and an unstable equilibrium (see Figure 9 ) . Obviously 
n = s + u. A secondary equilibrium class is an isomorphism class of quasi-duals, 
where the isomorphism is expected to preserve the colouring as well. 

While the derivation process above requires the existence of the gradient 
field of the height function, we would also like to consider some reasonable 
body surfaces with no gradient, e.g. a polyhedra. As the construction of the 
Morse-Smale complex is already extended for some non-smooth functions as 
well [TU] , we believe the definition of the quasi-dual and the secondary class of a 
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Figure 9: Columbus' algorithm: auxiliary coloured splitting Co defined only on 
the quasi-dual Pi of the Gomboc. 



polyhedron could be understood by intuition and needs no rigorous theoretical 
background. 

If a quadrangulation admits multiple non-isomorphic colourings, then its 
possible secondary classes form a partition of its map. It is easy to verify that a 
quadrangulation admits cither two non-isomorphic colourings, i.e. switching the 
colours results in another secondary class, or only one colouring. In the latter 
case we call the secondary class a self- dual. 

Assigning a quasi-dual to a convex body is a generalization of assigning a ra- 
dial graph to a polyhedral skeleton. Radial graphs are all quadrangulations, and 
every quadrangulation is the radial graph of some spherical surface graph [b"i|l6]. 
Moreover, assigning to a generic connected plane graph G its 2-coloured radial 
graph R(G) and colouring the images of the vertices of G "unstable" in R(G), 
and the images of the faces of G "stable" defines a bijection [TT]. For a fixed 
s, it, the bijection is between the set of secondary classes in the primary class 
{s,u} and the set of maps of generic connected plane graphs with u vertices, 
s faces and h = s + u — 2 edges. E.g. in the case of the minimal polyhedra, 
the quasi-dual is actually the coloured radial graph of the skeleton, thus the 
Poincare-Hopf theorem can be replaced by Euler's formula. 

The steps of Columbus ' algorithm are defined and their geometric feasibility 
were proved in |20) . where feasibility means they transform any convex body 
to another convex one. Observing the change of the quasi-dual of the under- 
lying body, the steps of Columbus' algorithm are combinatorial operators on 
quasi-duals we call the coloured splittings. A coloured splitting of a quasi-dual 
consists of a vertex splitting of the underlying quadrangulation followed by the 
proper colouring of the introduced vertex. We add that the coloured splittings 
corresponding to the original steps of Columbus' algorithm are all monotone. 
So each coloured splitting corresponds to two dual versions of a splitting of the 
quadrangulation: to one adding a new stable vertex, and to one adding a new 
unstable vertex. 

For technicality, we add an auxiliary coloured splitting Co applicable only 
on the quasi-dual P\ of the Gomboc, which is not literally a vertex splitting 
of a quadrangulation. Applying Cq either results in the quasi-dual in the class 



{1,2}, or in the quasi-dual in the class {2,1}, see Figure 9 So the definition 
of coloured splitting is relaxed to allow Co as well, considered a 1-splitting (i.e. 
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Figure 10: Hierarchy of secondary classes. Notation: X + : all secondary classes, 
Q + : generated from Pi, T: small secondary classes s.t. s + u < 8, X: irreducible 
ancestors, M: minimal polyhedra, A4 + : generated from minimal polyhedra. 



4.2 Generating secondary classes by monotone coloured 
splittings 

While the monotone coloured splittings generate every primary class from Pi, 
they do not generate every secondary class, so they admit a nontrivial family 
of irreducible ancestors. We say a secondary class is irreducible if it cannot 
be created with a monotone coloured splitting. In other words its underlying 
quadrangulation with n > 3 is irreducible, or n — 2. So the secondary classes 



in {1,2}, {2,1} are not considered irreducible (see Figure 9), but Pi in {1,1} 



is. By [Proposition 1| we can characterise the irreducible secondary classes for 
n > 3 as bodies with a quasi-dual of minimum degree 3. 



The results of this subsection are outlined in Figure 10 which we reveal 
gradually. Let Q denote the starting set consisting of the secondary class of 
the Gomboc, i.e. Q — {Pi}, and Q + the family generated from Pi by monotone 
coloured splittings. Throughout this section, for any starting set of secondary 
classes X, let X + denote the family of secondary classes generated from X 
by the monotone coloured splittings. So let X denote the family of irreducible 
secondary classes, then X + denotes the family of all possible secondary classes. 

The first part of |Corollary 1| states that X contains secondary classes besides 
Pi, which clearly follows from the fact that there is an infinite set of non- 
isomorphic irreducible quadrangulations. For the second part, let J- denote the 
family of secondary classes in the primary classes {s, u} such that s + u < 8. 



The second part states that T C G + , which follows from Theorem 2 with the 
help of the auxiliary splitting Cq. 



Note that Q + Hi contains only Pi by definition. Corollary 2 gives a geomet- 
ric characterisation of a special family within X called the minimal polyhedra. 
Recall that a polyhedron is a minimal polyhedron if its every face contains one 
stable and its every vertex is an unstable equilibrium, and let Ai denote the 
family of their secondary classes. E.g. the Platonic solids or a right prism are 
all minimal polyhedra, however, an oblique prism may not be a minimal poly- 
hedron: there may be a face which does not contain a stable point. |Corollary 2| 



states that M. c X (see Figure 10 1. As there are only three irreducible secondary 
classes until s + u < 10 according to our data set, we can safely say that every 
secondary class such that s + u < 10 are generated either from the Gomboc or 
from a minimal polyhedron, i.e. {S € X + : s + u < 10} C Q + U A4 + . 
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Figure 11: Irreducible subclass with parallel edges created by reflecting a simple 
one. 

Nevertheless, there are other irreducible secondary classes as well, e.g. |Fig- 



ure 11 shows a method to enumerate some of them with parallel edges. Clearly 



a minimal polyhedron cannot have parallel edges. 

Finally, [Theo rem 4| states that the irreducible ancestor of a secondary class 
is unique, because the quadrangulations of the irreducible secondary classes and 
the irreducible quadrangulations coincide for n > 3. Consequently, e.g. Q + and 
A4 + are disjoint. Theorem 4 also results in an algorithm to determine the ir- 



reducible ancestor of a body: we simply need to apply monotone contractions 
iteratively to its quasi-dual, until we get the irreducible ancestor. This algo- 
rithm is trivial because we can choose the contractions arbitrarily: the resulting 
irreducible ancestor is unique hence is independent of the actual choices, as any 
contraction sequence would result in the same ancestor. 



4.3 Generating primary classes by restricted splittings 

The restricted coloured splittings Sx t 2 (i.e. the monotone ones) can generate 
every primary class from the Gomboc [20j . but not every secondary class by 
|Corollary 1 We consider further restrictions in this subsection summarized 
by Corollary 3 In detail, the coloured splitting S±i generates every primary 



class from Pi, and the coloured splitting generates every primary class 
from the starting set of the secondary classes of the singleton primary classes 



section 5 



{{1, 1}, {2, 1}, {3, 1}, {1, 2}, {1, 3}}. Note that according to |Table 1 
these are the only singleton primary classes. 

Proof of Corollary 3. First we consider restricting the monotone coloured split- 
tings to S^i. As 1-splitting is always applicable at any part of the graph 
if n > 2, any primary class can be reached with coloured 1-splittings from 
{1, 2} or from {2, 1}. So the primary classes are generated from the starting set 
{{1, 1}, {1, 2}, {2, 1}}, and this set can be generated from Pi by Co- 

Now we consider restricting the monotone coloured splittings to 52,2- The 
quadrangulations of the classes {2, 1}, {3, 1}, {1, 2}, {1, 3} are 2-irreducible (see 



Figure 5), so they has to be in the starting set. The 2-splitting is always appli- 



cable around a vertex v with d(v ) > 2 introducing a vertex w of the same colour 
as v. It is clear that if n > 3, both independent sets of a bipartite quadrangu- 
lation contains vertices of degree at least 2. As C4 in {2,2} is 2-contractible, 
the 2-splittings can be arbitrarily combined to generate some secondary class in 
any primary class {s, u} such that s + u > 4. □ 
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In the next section we show statistics on the number of 1-irreducible and 
2-irreducible secondary classes for a limited size in Table 3 



5 Computational results 

This section presents some statistics on the data set attained by the computer 
program. As plantri supports dividing the computation into independent 
parts, we could perform this computation in parallel in a grid infrastructure 
using the Saleve framework |5|. Some of the numbers have been already pre- 
sented in 131. 



Table 1 shows the cardinalities of the classes {s,u} with s + u < 10. The 
table is symmetric because a subclass such that s ^ u is clearly not self-dual. 
These numbers were already published as the cardinalities of the unrooted and 
unsensed maps by Wormald [24] , Walsh [2TJ [23] . 



Table 2 shows the number of multiquadrangulations q(n), which is the first 
publication of these numbers, up to our best knowledge. To compute q(n) from 
the same program-generated data set, observe a relation between the number of 
secondary classes in the class {s,u}, the number of self-dual secondary classes 
and the number of quadrangulations of size n — s + u and, denoted respectively 
by e(s,u), e SD {s,u), q(n): 

n-1 

2g(n) - e SD (n/2, n/2) = £ e(s, n-s). (1) 

|Table 3] shows the number of ancestor secondary classes with respect to 
splittings with different restriction criteria: 52,2, <Si,i and Sj.,2. 
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s = 1 
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Table 1: Cardinalities e(s,u) of the equilibrium classes. 
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30 
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52 
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= 7 


124 




248 


n 


= 8 


733 


50 


1416 


n 


= 9 


4586 




9172 


n 


= 10 


33373 


380 


66366 



Table 2: The number of multiquadrangulations (q), self-dual secondary classes 
(csd) and secondary classes (J2 e )- 







1-2-contractible 


2-irrcduciblc 


1-irrcduciblc 


irreducible 
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= 4 





3 
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= 5 


6 


2 


6 





n 


= 6 


32 


4 


16 





n 


= 7 


172 


10 


66 





n 


= 8 


1071 


33 


311 


1 


n 


= 9 


7370 


114 


1688 





n 


= 10 


55766 


474 


10125 


1 



Table 3: Ancestor secondary classes with respect to S^i (third column), 62,2 
(second column), S\ t 2 (fourth column) and non-ancestors (first column). 
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